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We suggest a structure for the vacuum comprised of a network of tightly knotted/linked flux tubes 
formed in a QCD-like cosmological phase transition and show that such a network can drive cos¬ 
mological inflation. As the network can be topologically stable only in three space dimensions, this 
scenario provides a dynamical explanation for the existence of exactly three large spatial dimensions 
in our Universe. 


Introduction: Although the question of why our Uni¬ 
verse has exactly three (large) spatial dimensions is one 
of the most profound puzzles in cosmology - especially in 
view of quantum gravity scenarios such as string theory 
which assume 9 or 10 space dimensions at the funda¬ 
mental level - it is actually only occasionally addressed 
in the literature [l]]. In this paper we propose a topo¬ 
logical explanation for the dimensionality of space-time 
based on the idea that inflation is driven by a tightly 
knotted network of flux tubes generated in a cosmolog¬ 
ical phase transition and the fact that one-dimensional 
knots are topologically stable only in exactly three space 
dimensions (’’Knotty Inflation”). 

The main idea of the model is connected to the fact 
that, in non-Abelian gauge theories analogous to QCD, 
chromoelectric flux tends to become confined and the re¬ 
sulting tube-like structures can be treated as effective 
one-dimensional objects or strings. The formation of 
flux tubes is behind the basic models of hadronization in 
QCD, with flux strings connecting quark-antiquark pairs 
giving rise to a linear potential that confines them into 
mesons and prevents the existence of free colored charges. 
A similar description with more complicated flux tube 
shapes also describes the formation of baryons and glue- 
balls, the latter denoting bound states of pure flux. 

During the phase transition from a generic quark-gluon 
plasma to a hadron gas, a large number of flux tubes may 
fill up the whole Universe and form an intricate network 
of both open and closed flux tubes, with some similarities 
to what is known as the ‘spaghetti vacuum’ (2;,], poten¬ 
tially with a large linking and/or crossings number den¬ 
sity provided the flux tube density is sufficiently large. 
A related idea has been suggested recently in Q. It has 
been shown that such a network is stabilized by a topo¬ 
logical conservation law @]. 

This network will then tend to relax into a tight con¬ 
figuration that fills the whole of space in order that each 
tube minimizes its energy, similarly to the description of 
glueballs in terms of knotted/linked flux tubes proposed 
in Q. The network then decouples from the Hubble flow. 

The energy stored in this tight network will thus pro¬ 


vide an effective cosmological constant which is homoge¬ 
neous on scales large compared to the typical flux tube 
width and on average isotropic if the flux tubes are ran¬ 
domly oriented. This effective cosmological constant will 
thus trigger a period of inflation in the early universe that 
lasts until the network decays through reconnection, tube 
breaking and other quantum effects. Network breaking 
will proceed until a gas of hadrons and radiation is formed 
and the standard cosmological evolution begins. 

In higher-dimensional theories, such as string theory, 
gauge theories are typically confined to hypersurfaces of 
p spatial dimensions known as Dp-branes, to which open 
strings with gauge, fermionic and scalar degrees of free¬ 
dom are attached. These fundamental strings behave 
like point-like objects on energy scales smaller than the 
fundamental string scale, M s , and so, for a confinement 
scale A < M s , we nevertheless expect the formation of 
tube-like structures connecting color charges in higher¬ 
dimensional gauge theories. The tube cross sections will 
have a higher-dimensional geometry, in the simplest case 
with S p ~ 2 topology, but will in any case behave effec¬ 
tively as one-dimensional objects. 

The main difference between the dynamics of confine¬ 
ment in 3 + 1 and higher dimensional gauge theories lies 
in the fact that flux tube knots and links will only be 
topologically stable (or metastable) in three spatial di¬ 
mensions, being otherwise able to unknot/unlink in the 
extra-dimensions. Given that this is an essential feature 
for the formation of a tight network of flux that decouples 
from the Hubble flow and provides an effective cosmolog¬ 
ical constant, confinement can only lead to inflation for 
gauge theories in D3-branes. 

Hence, considering a generic universe in a higher¬ 
dimensional space with gauge theories living on branes 
of different dimensionalities, only D3-branes will inflate 
with the formation of a tight flux tube network during 
a cosmological phase transition. This can thus explain 
why our universe has only three large spatial dimensions 
in the absence of other sources of vacuum energy. 

It is the purpose of this Letter to outline the main fea¬ 
tures of this mechanism and discuss the necessary prop- 
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erties that a flux tube network must exhibit in order to 
yield a successful inflationary model in three dimensions. 
While leaving a detailed modeling of the network’s evo¬ 
lution and observational predictions of this scenario for a 
future work, we will describe the most important physical 
processes that drive the network’s formation, inflationary 
dynamics and subsequent decay, giving quantitative es¬ 
timates for the relevant time and energy scales. 

Flux tubes in gauge theories: Understanding confine¬ 
ment in QCD and gauge theories in general is one of 
the most important problems in particle physics. Due 
to the strong coupling and non-perturbative nature of 
confinement, this is an intrinsically hard problem that 
can currently only be studied accurately in the context 
of numerical simulations on the lattice. It is nevertheless 
widely accepted that the confinement of chromoelectric 
charges, such as quarks and anti-quarks in QCD, is asso¬ 
ciated with a squeezing of chromoelectric flux lines into 
tube-like structures connecting the charges. These flux 
tubes then give rise to a potential energy that grows lin¬ 
early with the distance between pairs of quarks and anti¬ 
quarks, therefore leading to their confinement in mesons, 
baryons and potentially other hadronic structures. Flux 
tubes can be described e.g. in terms of Abrikosov-Nielsen- 
Olesen vortices [6, M] in dual superconductor models of 
confinement (see [ 8 j for a review). 

Reference [fjj pointed out that closed flux tubes may 
have non-trivial topologies, including knotted tubes and 
links between distinct closed flux tubes. The degree of 
“knottedness” of each configuration is then associated 
with a conserved topological charge, analogous e.g. to 
the (Abelian) magnetic hclicity that characterizes con¬ 
ducting fluids in magnetohydrodynamics Q. Any such 
topologically non-trivial flux tube configuration will nec¬ 
essarily relax into an equilibrium state by minimizing its 
length and therefore its energy. These equilibrium states 
then correspond to the tightest knotted or linked config¬ 
urations with a given topological charge and number of 
flux quanta. This thus motivates us to consider the de¬ 
scription employed in jf|, where the chromoelectric field 
F(a (where i = 1,2,3 in three spatial dimensions) is con¬ 
fined in knotted/linked tube-like structures that carry 
one quantum of chromoelectric flux. The relevant (static) 
Lagrangian density is given by: 


£ = ±tr[F 0i F oi ] + trA 


&E 

( to 2 ) 



(1) 


up to the addition of a constant energy density as in the 
MIT bag model jToi], where the last term enforces flux 
conservation across the tube sections through a Lagrange 
multiplier A, with n* denoting the unit vector normal to 
the tube’s cross section. We expect the flux tube radius a 
to be given by the confinement scale A, such that a ~ A -1 
in natural units, which we will employ in our discussion 
henceforth. This yields the equations of motion for the 


gauge field: 

D°(F 0i - A m) = 0 , D l (F 0i - A m) = 0 , (2) 


the solution of which is given by a constant chromoelec¬ 
tric field 

Foi = — (3) 

TO 

which vanishes outside the flux tube. The energy density 
inside each flux tube is then given by: 


1 tr$|. A 4 

2 (to 2 ) 2 27t 2 ’ 


( 4 ) 


up to 0 ( 1 ) factors that will not affect our discussion. 

Network formation: Let us then hypothesize that there 
exists a confining gauge theory, analogous to QCD but 
with a generic high-energy confinement scale A, and 
which lives on a three-dimensional brane within a higher¬ 
dimensional compact space. For simplicity, we will de¬ 
note the particles charged under this gauge group as 
quarks and gluons, although one must bear in mind that 
these are not the known QCD fields but rather novel de¬ 
grees of freedom that must decouple from the low-energy 
effective theory. 

In the early Universe, for temperatures above a cer¬ 
tain critical value T c ~ A, these quarks and gluons are 
essentially free, forming a plasma that we will assume is 
close to thermal equilibrium and dominates the energy 
balance in the Universe. This plasma has an energy den¬ 
sity Pr(T) = (7r 2 /30)g*T 4 , where g * denotes the num¬ 
ber of relativistic degrees of freedom, which includes at 
least the massless gluons and the light standard model 
fields of the theory. Due to expansion, the temperature 
of the plasma will decrease until it becomes lower than 
the critical value, at which point all the free gluons be¬ 
come confined in both open and closed flux tubes, the 
former connecting quark-antiquark pairs. Given that 

Pr(T c ) ~ 3 * A 4 - p E , (5) 


we expect that, at the confining phase transition, a large 
number of flux tubes is formed within each Hubble vol¬ 
ume. Such a large density will naturally lead to a large 
number of crossings between the flux tubes and, hence, 
to a large number of knots and links of different configu¬ 
rations. This is analogous to the behavior observed when 
a string is tumbled inside a box of fixed volume, where 
the probability of knotting increases very quickly above 
a certain critical string length and which effectively cor¬ 
responds to a critical density [11] ■ The flux tube network 
will then be endowed with a non-trivial topology, which 
can be metastable in three spatial dimensions and that, 
as we describe below, can be sufficiently long-lived to 
drive a period of inflation in the early Universe. 

Inflation: After the phase transition, Hubble expan¬ 
sion will tend to stretch the flux tubes, while maintaining 
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their radius (set by the strong non-perturbative dynam¬ 
ics), and increase the size of the gaps in between them. 
However, as argued above, the presence of knots and links 
makes the network behave differently from a system of 
isolated flux tubes, endowing it with a rigidity and mak¬ 
ing it try to relax into a tight equilibrium configuration. 
Flux tubes will then shrink and approach each other in 
order to maximize the fraction of the spatial volume that 
they can occupy given the topological constraints. In ad¬ 
dition, the number of knots and links in the network will 
decrease due to the decay processes that we describe be¬ 
low. These will make the network progressively less tight, 
with the tightest configuration occupying an increasingly 
smaller fraction of the spatial volume. 

We then expect a network to remain tight if: 

T~tight H T~decay 5 (b) 


i.e. if it relaxes into a tight configuration more quickly 
than expansion and knots/links are stable on the Hub¬ 
ble time scale. If the network is sufficiently dense at 
the phase transition, i.e. if the gaps between flux tubes 
are not much wider than the tube radius, we expect 
Ttight A -1 , and the network can in principle relax into 
a tight configuration within a Hubble time for A>il. In 
this case, the tube width is also much smaller than the 
Hubble radius, so that the network is homogeneous on 
near-horizon and super-horizon scales. It should also be, 
on average, isotropic if quarks and gluons are randomly 
distributed in the thermal plasma at the transition. 

Let us see that a tight network will behave as an effec¬ 
tive cosmological constant as long as the above time scale 
hierarchy is satisfied. The “tightness” of the network can 
be evaluated in terms of the fraction of the spatial vol¬ 
ume occupied by the flux tubes, /, which will in general 
depend on the space and time coordinates. On average, 
we expect a network that is homogeneous and isotropic 
on scales larger than the tube radius, such that / = fit) 
at leading order. The network thus has an average energy 
density p = Jpe, with the first law of thermodynamics 
yielding: 

d iP v ) = ^(PEV)dV + fp E dV = -pdv , (7) 

neglecting decay processes and consequent heat transfer. 
From this we conclude that the network has an effective 
pressure: 


P=~ 1 + 


1 dlog/N 
3 dN e ) 


( 8 ) 


where N e denotes the number of e-folds of expansion, and 
hence an equation of state parameter w = p/p ~ —1 for 
/// « 3 H. 

The network will initially relax into a tight configura¬ 
tion with dlog // dN e > 0, yielding a phantom-likc equa¬ 
tion of state w < — 1 [HJ . This phantom behavior is often 


considered in dynamical dark energy models, as well as 
in dark matter scenarios with bulk viscosity [13| , and 
it would be interesting to consider this behavior in the 
context e.g. of axionic strings 14]. 

The network will then subsequently remain close to a 
tight configuration, with a fixed or at most slowly vary¬ 
ing filling fraction /. The rigidity of the network, which 
results from the presence of a large density of knots and 
links, will thus decouple it from the Hubble flow and gen¬ 
erate an effective cosmological constant. 

As the number of knots and links decreases due to net¬ 
work decay, / will then decrease such that w > — 1, thus 
yielding a period of inflation with an equation of state 
analogous to that of a canonical slowly-rolling scalar field. 
Note that, in the absence of the topological charge that 
maintains the tightness of the network, the filling frac¬ 
tion would decrease quickly and no accelerated expansion 
could be obtained. 

The Friedmann equation yields, during inflation, 


H 2 ~ 


/A 4 

67T 2 Af2 


(9) 


which implies that A < H for a sub-Planckian confine¬ 
ment scale, with e.g a phase transition at the grand uni¬ 
fication scale, A ~ 10 16 GeV, yielding H/A ~ 10 -3 for 
/ < 1. This shows that the flux tube network will in gen¬ 
eral be on average homogeneous and isotropic on super¬ 
horizon scales and that relaxation may in principle occur 
in less than a Hubble time, as required above. 

Network decay and reheating'. The flux tube network is 
metastable in three spatial dimensions, since its degree of 
knotting/linking is reduced by string breaking and other 
reconnection events. 

Two flux tubes may unknot or unlink by quantum tun¬ 
neling. We can estimate this by treating the intersection 
between two flux tubes as a (non-relativistic) particle of 
mass m ~ p E a 3 ~ A/2ir 2 that tunnels through a po¬ 
tential barrier of length ~ a and height AV > A. The 
standard tunneling rate computation then yields a typi¬ 
cal lifetime: 


1 , . I 

r t H = —{aH) exp - 
P \ 7T 



( 10 ) 


where /3 is the typical flux tube velocity. The tunneling 
rate is thus suppressed for a sufficiently large potential 
barrier separating the flux tubes. Such a barrier has been 
obtained for SU(2) flux tubes within the dual supercon¬ 
ductor description of confinement [15J but, as far as we 
are aware, more generic computations have yet to be de¬ 
veloped. In bag models, the flux enclosed within a tube 
is supported by surface currents and one can expect the 
repulsion between currents in adjacent tube sections to 
effectively provide a ‘hard’ boundary for the flux tubes, 
making tunneling more difficult. 
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If tunneling is indeed suppressed, flux tubes will pri¬ 
marily break, and hence unknot or unlink, through the 
well-known Schwinger effect fl6| , where QQ pairs are pro¬ 
duced in the approximately constant chromoelectric field 
enclosed within each flux tube. For Nf quark flavors of 
mass toq, using that the string tension k = ira 2 pE ~ 
A 2 /(27t), we obtain for the characteristic lifetime of a 
flux tube of length l: 


r b H 


8tt 4 


N 


f 


(aH) 2 (lH) 4 exp 27 t 


ln Q 

A 2 


( 11 ) 


This implies that the lifetime for string breaking can be 
parametrically large if there are no light quarks in the 
spectrum of the gauge theory, toq A, even for a <C 
H~ x and flux tubes as long as the Hubble radius. Thus, 
the confining gauge theory responsible for a sufficiently 
long period of inflation in the proposed scenario has to 
be distinct from QCD, since in the latter case knots and 
links between flux tubes will quickly decay through the 
production of up and down quark-antiquark pairs. 

Given the exponential suppression of the string break¬ 
ing and tunneling rates, it is thus not difficult to envisage 
scenarios where the network remains tight for the 50-60 
e-folds of accelerated expansion required to explain the 
present flatness and homogeneity of the Universe. 

String breaking will eventually lead to a system of un¬ 
knotted/unlinked flux tubes, and thus a gas of hadrons. 
These may be unstable and decay into radiation, includ¬ 
ing in principle the Standard Model particles. The de¬ 
tails of this process depend, of course, on how the con¬ 
fining gauge theory at a high energy scale A is related to 
the low-energy physics, but it is clear that the proposed 
model naturally includes a mechanism for reheating the 
Universe after inflation. It is also possible for a significant 
amount of radiation to be produced during the inflation¬ 
ary period, if e.g. small loops unlink from the network 
and decay while the network is still tight, thus poten¬ 
tially leading to a warm inflation model 17]. 

Cosmological perturbations: The flux tube network is, 
on average, homogeneous on scales larger than the tube 
radius, which as shown above can be parametrically be¬ 
low the Hubble radius. Nevertheless, fluctuations in 
the energy density of the flux tube network may arise 
on super-horizon scales and become imprinted on the 
curvature of space-time, thus seeding the temperature 
anisotropies in the Cosmic Microwave Background and 
the Large Scale Structure of the present Universe. 

Curvature perturbations in a perfect fluid driving a 
period of accelerated expansion propagate with an imag¬ 
inary sound speed c 2 = dp/dp = w < —1/3, which would 
make the system unstable to small fluctuations as op¬ 
posed to the canonical scalar field models of inflation. 
However, the flux tube network has a non-vanishing rigid¬ 
ity, as a result of the presence of numerous knots and 
links. The network thus behaves more like an elastic solid 


than a perfect fluid, similarly to the networks of topolog¬ 
ical defects considered in 18, M- The anisotropic stress 
inherent to elastic solids allows for the propagation of 
both longitudinal and transverse waves and the associ¬ 
ated sound speed c 2 = w + (4/3)p/(p +p) can be real for 
a sufficiently large shear modulus p [211 . 

Inflationary models with elastic solids may lead to a 
nearly scale-invariant spectrum of both curvature and 
tensor perturbations, with interesting differences from 
scalar field models [20(. Firstly, despite the absence 
of non-adiabatic modes, perturbations evolve on super¬ 
horizon scales due to the presence of anisotropic stress, 
although maintaining the relevant scaling between differ¬ 
ent super-horizon modes. The overall amplitude of the 
spectrum thus depends on the details of reheating, in our 
case of the decay of the flux tube network. Anisotropic 
stress will also act as a source for tensor perturbations, 
potentially modifying the primordial gravitational waves 
spectrum with respect to canonical models. 

Although concrete observational predictions require a 
detailed modeling of the properties and dynamics of the 
flux tube network, which is beyond the scope of this Let¬ 
ter, it is worth mentioning that a red-tilted curvature 
spectrum can only be obtained in elastic solid models of 
inflation for a slowly varying equation of state. We expect 
this to be the case in our model since a knotted network 
will be close, but not exactly in, a tight configuration, 
as a result of the opposing effects of Hubble expansion, 
relaxation and unknotting events. 

Conclusion: A knotted/linked flux tube network 
formed in a QCD-like phase transition can provide a 
natural source of inflation. Furthermore, this may ex¬ 
plain why we live in three large spatial dimensions, since 
knotted/linked one-dimensional tubes are topologically 
unstable in higher-dimensional space-times. We realize 
that the picture developed in this paper may also be ap¬ 
plied to a model of dark energy, which would eliminate 
the need for an ultra-light scalar field. 

The volume filling fraction / is not a fully dynami¬ 
cal variable because it does not enter the field theory 
Lagrangian in the above calculations. One can easily 
(and most likely without any unpleasant consequences) 
remedy this drawback by considering a fully dynami¬ 
cal model. As a concrete, simple and non-trivial ex¬ 
ample of such, we can take the abelian Nielsen-Olesen 
model, where the role of the fraction / is played by time- 
dependent characteristics of spatial extents of the scalar 
and gauge fields. It is plausible that main dynamical 
characteristics of flux tube knotting and the evolution 
of a knotted network in terms of tightening, Hubble ex¬ 
pansion and decay in this abelian model will not change 
significantly when a generic non-abelian gauge theory is 
used instead. Although exact solutions in such models 
are unavailable, appropriate approximations should be 
enough to establish quantitative features to support the 
rest of our calculations. We will present the detailed anal- 
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ysis based on the abelian Nielsen-Olesen model elsewhere 

Ill- 
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